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2-D Nonlinear Theory of the Free 
Electron Laser Amplifier for an 
Electron Beam with Finite 
Axial and Transverse Dimensions 

Many current experiments of the free electron laser (FEL), ultilize electron beams from a 
millimeter to a few centimeters in pulse length. The short pulse length is typical of high energy 
accelerators such as RF Linacs and microtrons. The finite length effect of the electron beam on 
the radiation was found to be important in the Stanford oscillator experiment.* 11 Currently, 
many experiments in the amplifying mode are being conducted with the short electron beam 
pulses. We have an analytical expression for the gain pulse of the radiation field, applicable to 
these experiments. 

One-dimensional analysis of the radiation field for electron beams of finite length have 
been numerically simulated on computers.* 241 The effect of the finite transverse dimensions 
was either not included, or incorporated through filling factors. On the other hand, previous 
three-dimensions self-consistent formulation* 51 of the radiation field for a semi-infinitely long 
electron beam in the amplifying configuration has resulted in a number of interesting effects 
not obtainable by the 1-D formulation. Numerical effort to find the growth of the 3-D radia¬ 
tion field on the finite length electron beam began only recently.* 61 . 

In this paper, we will present a fully 2-D, self-consistent, non-linear, analytical analysis of 
the FEL process in the amplifier mode of operation treating the finite length and transverse 
effects associated with both the electron beam and the radiation beam. Our formulism also 
includes various efficiency enhancement schemes: (i) contouring in the longitudinal direction 
the amplitude and/or the wavelength of the magnetic wiggler field, and (ii) applying an external 
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D.C. electric field. Analytical results in the amplifying configuration are obtained in the low 
gain, trapped particle regime. 

The schematic of the configuration is shown in Fig. 1. The generalized vector potentials 
of the right-handed, heliccl, static magnetic wiggler field and the electromagnetic radiation held 
are 

A„ (z) « /f„ (z) [cos ( f k w (z') dz) e x + sin ( f k w (z')dz‘)eA M) 

•/0 Jo 

\r (x. y. z, t ) - Ao (x. y, z. I) [cos (— z - m i -+ p (x. y. z. i)) e x 

c 

— sin (— z - to / + <t> (x, y, z, /)) e.l (2) 

c 

where A w and k w are all slowly varying amplitude and wave number of the wiggler field and A R 
and >p are slowly varying amplitude and phase of the electromagnetic radiation field following 
the electron pulse. We also include an external DC electric field, E oc (z) - -d<t>oc(z)/dz e 2 
for the purpose of efficiency enhancement. 


In this analysis we will not consider the gradient in the wiggler field. This is a good 
approximation if k w r b « 1, where r b is the radius of the electron beam. If the FEL is operat¬ 
ing in a trapped particle mode, we also require 15 ’ r b < (y z0 O -1 (8 V2 y z a/|8 a ) l/2 (A R IA„) u * y 
where y z0 - (1 - v 2 0 /c 2 )~ in , /3 a - Id Aj(y 0 m 0 c 2 ), y 0 - y z0 ya., ya “ 

(1 + |e I 2 Aj(0)/(mo c 4 )) l/J , and v z0 is the axial velocity at z ** 0. 


The electron motion can be described in terms of their phase <ti in the ponderomotive 


1 d 2 >fi _ _1 _d 2 f(z,t) + K (z) 

\ u> 1 ( lei ] 2 9><»(z) 

c 2 dt 2 c 2 dt 1 9z 

2 f y 2 l "» 0 c 2 ) 


(O 1 

Id 

ddtQQd) 

, 2k w (z) 

Id 

c yy z 2 

m 0 c 2 

dz 

- y 2 r 

m 0 c 2 


A w (z) A k sin 1 1 


(3) 


'•Wo-fo-') 


- C Z V’ V -IT 

where if»(x 0 , yo- (o> t) - J 0 (/c w (z') + <o/c)dz' + a >t + 'fix, y, , z, t) is the phase for 

the electron, which was at (x 0 , .yp, (o) at / - 0, y - y t y L . y z - (1 - v?/c 2 )~' n , y x - 
(1 + ld J /4»(z)/(m^c 4 )) l/2 , v z - Idi/dt- d<e!dt + cel/U^fz) + w/cl is the axial velocity, z — 
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*o + /„ Vj (x 0 , > 0 , lo, t')dt' is the axial position of the electron and f 0 is the axial position of 
the electron relative to the center of the electron beam at t - 0. 

The wave equation for the radiation field is (V 2 - c~ 2 d 2 /dt 2 )A R - -4flrc -, J , where 

j » — H-° J* rf| 0 J" ^ dxo f _ ** (*o. >o)* Cfo) 

m 0 

-y _l Pi8(x - x 0 )8(y - y 0 )8(z - z) (4) 

is the current, (x 0 , yo) are the particle's transverse positions at t — 0, 0 (x 0 , > 0 ) is the 
transverse current density profile, > s the macroscopic electron pulse shape, n 0 is the peak 
(€ I 

current density, and P L — — A w is the transverse momentum. 

c 


We can rewrite the radiation field as A* — a R (x. y, z, t) exp [i(uiz/c — ojf)]e+ + c.c., 
where a, — A R exp(i<p) is the complex amplitude of the radiation field, and e ± ” (e x ± ie y )/2 
is a new coordinate system. The wave equation assumes the form 


+ + f + if 

dx 2 by 2 c dz c dt 


a R (x, y, z, t) 


~ J*_. d t° J-oo <**<>/_. Cfo. j'o) A (f o) 


/4 w (z) 


exp|-ijj^ (A: w (z') + <u/ c)dz'— wfJJs(x — x 0 )6(y- y 0 )6(z- *)• (5) 

In obtaining (5), we have used the fact that l(3 2 /3z 2 - c~~ 2 b 2 !bt 2 )a K I « 2a>/e|(3/3z + 
c~'d/dt)a R \. 


The solution for a R can be separated into the input radiation field a ( „, and the excited 
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The integral in time of Eq. (6) can be evaluated by changing the argument of the delta func¬ 
tion. 


“ Jo dt ' f_ m d i o dx 0 f_" dyof(x 0 . y Q , f 0 . x. y, z, /') —-— 

dt' 


^ /(-yQ. J>0. gp, X, y. z. T 0 ) 

° C- v z (x 0 . y 0 , fo, r 0 ) 


where r 0 (x 0 . yo, io, t) is the retarded time associated with the electron, which originated at 
(•*<>• yo- it) at t ” 0. The retarded time satisfies the equation 


to + z f (r 0 > + c(t - r 0 ) = £ + z c (t) (8) 

where z c (t) — \ sc (t')dt' is the macroscopic location of the center of the electron beam at 

time /, Vj c (t) — <o/(k w (z c {t)) + <u/c) is the macroscopic velocity of the electron pulse, f is the 
position of the electron relative to the center of the electron beam at time t. 


The complex radiation amplitude in Eq. (8) can be evaluated if we make the following 
simplifying assumptions. For experimental parameters of interest, we can assume that the 
bunching mechanism does not alter the macroscopic electron pulse shape, hence, it travels 
undistorted through the interaction region. We will assume that the electron beam has an axi¬ 
ally symmetric Gaussian profile in the transverse direction, i.e., 0 (x 0 . yo) “ 
exp[-(x 0 2 +yo 2 )/'’o J J- Furthermore, we will assume that the waist of the input radiation field r 0 
is much larger than the radius of the electron beam r b , such that i is approximately a function 
of ( 0 and / only. The excited radiation field takes the form 

, A w (j + z c (t ) - c(t- r 0 )) (1 + v,/c)y/ 

y cit — r 0 ) — iz b (9) 

e*P i c(f _ _ /z J | exp (-/(»fr(^o, r 0 ) + f, t 0 ))1 
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where z b - r b <a/2c is the Rayleigh length associated with the electron beam radius. Equations 
(3) and (9) describe self-consistently a general, nonlinear, 2-D, FEL amplifier with a macros¬ 
copic pulse shape A(f 0 ). 

For the purpose of illustrating the finite length pulse effects in an FEL amplifier operating 
in the low gain limit, i.e., |a in | » |a„|, we take the electron beam profile to be uniform, i.e., 
h(g 0 ) * 1 for |f 0 | < L b /2 and h(to) = 0 for |f 0 l > L*/2, where L b is the length of the elec¬ 
tron pulse. We also make the constant phase, resonant particle approximation. In this approxi¬ 
mation all particles are assumed to have the same constant phase, i R . The electron beam in 
this approximation consists of a pulse train of macro particles separated in distance by 2nv 0z /b». 
Furthermore, we will limit ourselves at this point to a constant parameter wiggler and consider 
only an external DC electric potential. The amplitude and phase of the total field are 

A R (r, t . l) =» A,„ - a$A w \1, cos' ~l> R + /, sin Ji R ) (10a) 

<f(r.t. t) = -ao<4»I/, cos i ]> R - l r sin<i»*] (10b) 

where A,„ — |a in |, I r — Re (/); /, — Im (/), / * £;|—y-«/J - £,|—£, is the 

exponential integral function, 

<?i- (-iz»)(2-y 2 R (f - to.,) - iz b )'\ 

Q u “ (~i 2 b )( 2 y} R (t - to.J -*»r\ 

to.u ” f (for ( < LJ2) and t “ L b /2 (for t ^ L b /2) is the upper limit of the integration, to.i 

— t -(c-v K )t (for t ~(c- y K )t > -L b /2) and to.i “ ~L b /2 (for t -(c- v w )r < -L b /2) is 
the lower limit of the integration, and y zR is the resonant gamma associated with the axial 
motion. 

A more realistic electron beam profile (i(f 0 ) *■ 1 -(2t^L b ) 2 (for to ^ L*/2) and h(t o) 

- 0 (for If 0 | > LJ 2) can also be integrated. The result is not given here, because the more 
complicated expressions would obstruct the initial understanding of the physical process of the 
pulse propagation. 
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As an example of a 10.6 /um FEL utilizing a C0 2 laser as an input field, we choose an 
electron beam of energy 25 MeV (y 0 — 50), current of / - 5 A and radius (Gaussian profile) 
of r b — 0.5 mm and pulse length L b — 3 mm. Such a beam has a peak density on axis of n 0 
— 1.3 x 10" cm -3 ( <u b - 2.0 x 10 10 sec _l ). The constant parameter wiggler has a magnitude 
of B w — 5.0 kG and wavelength of /* = 2.8 cm which gives A w - 2.2 x 10 3 statvolts. The 
wiggler velocity is v 0 i “ 2.6 x 10 _J c which gives “ 1.35 and y t * 37. The input C0 2 
power density is taken to the P,„ *= 4 x 10* W/cm 2 which gives A ln — 0.30 statvolts. Our illus¬ 
tration assumes resonant macro particle approximation and an applied D.C. electric potential 
such that sin ^ R = 0.6. 

The schematics of the gain 

G(r.|. ()- (Axir.t. t) - AJ!A in 

are shown in Figs. 2 and 3. The slashed bars in the (z. t ) plot of Fig. 2 denote the locations of 
the electron beams at t x — 1 m/c and t — 2 m/c, which c is the speed of light. The solid lines 
in the (z, f) plot are the light lines. The gain pulse on axis are plotted at times f, and r 2 . We 
see that the excited radiation pulse grows and spreads beyond the electron beam pulse. The 
transverse variation of the gain at £ — 0 for various times are plotted in Fig. 3. The decrease 
of radiation field far from the axis is due to refraction toward the center of the beam. 

We have obtained a general expression for the growth of the 2-D, stimulated radiation 
pulse on an electron beam of finite axial and transverse dimensions in an FEL amplifier. We 
included diffraction as well as refraction. In the axially symmetric, low gain, resonant macro 
particle limit, we have an analytical expression for the radiation gain. The formalism presented 
here can be modified to study the radiation build up and ” laser lethargy” in the FEL oscillator. 
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